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Abstract. The fluid dynamic limit of the Boltzmann equation leading to 
the Euler equations for an incompressible fluid with constant density in the 
presence of material boundaries shares some important features with the better 
known inviscid limit of the Navier-Stokes equations. The present paper slightly 
extends recent results from [C. Bardos, F. Golse, L. Paillard, Comm. Math. 
Sci., 10 (2012), 159-190] to the case of boundary conditions for the Boltzmann 
equation more general than Maxwell's accomodation condition. 



1. The Inviscid Limit of the Navier-Stokes Equations with Dirichlet 

Condition 

Consider the flow of an incompressible fluid with constant density confined in 
a domain £7 C R^ (in practice, = 2 or 3), with smooth boundary dQ. The 
dimensionless form of the Navier-Stokes equations governing the velocity field of 
the fluid Uf: = Ue{t,x) G R^ is 

' divx tie = , 

dtUg + divxiuf^ (8) Ue) + '^xPe — eAa;Mg , X e , i > , 



(1) 



where e = Re^ is the reciprocal Reynolds number of the flow. The Dirichlet 
boundary condition Ue | = means that the fluid does not slip on the boundary 
5r2 of the domain. 

An outstanding problem in fluid dynamics is to understand the behavior of in 
the vanishing e limit, and especially to decide whether, in that limit, converges 
to the solution of the Euler equations 

( div,r u = , 

X e fi, t> , 



(2) 



dtU + divx (u (8) u) H- VxP — , 



0, 



it=o 



where n is the unit outward normal field on dfl. 
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In the case of the Cauchy problem set in cither the whole Euclidian space R or 
the periodic box T''^, i.e. in the absence of material boundaries, the Navier-Stokes 
solutions converge to the solution of the Euler equations as long as the latter remains 
smooth — see for instance Propositions 4.1 and 4.2 in [21]. We recall that, if N = 2 
and e C^'" (i.e. either C'^'"(R^) or C'='"(T^)) for some fc > 1 and a € (0, 1), 
there exists a unique global solution u G C(R+;C'^'") of the Euler equations — 
see for instance Theorem 4.1 in [21] — while, if A'' = 3 and if u*" € C^'", there 
exists a unique local solution u € C{[0,T);C^'°') for some T € (0, +oo], and it is 
unknown at the time of this writing whether T = +oo. Thus in the absence of 
material boundaries, the inviscid limit of the Navier-Stokes equations is described 
by the Euler equations, globally in time if A'' = 2, and maybe only locally in time 
if iV = 3. 

The situation is completely different in a domain with material boundaries. Even 
when (2) has a smooth solution, this solution might fail to attract u,. in the limit 
as e ^ 0. The reason for such a behavior is of course that the Dirichlet condition 
satisfied by the Navier-Stokes solutions Ue for each e > is overdetermined for the 
Euler equations. Indeed, in general, the tangential component of the solution u of 
the Euler equations does not vanish on d^l. 

Thus one may seek to match the Euler solution to the Dirichlet condition on the 
boundary with a viscous boundary layer of thickness -^/e as proposed by Prandtl 
— see for instance [29] . However it may not always be possible to do so. Consider 
for instance the flow of a viscous fluid past a cylinder or a sphere. In other words, 
assume that is the complement in R^ of a circular cylinder — or the complement 
in R^ of a ball, and that the velocity field is constant at infinity on one side of the 
immersed body. Although there is no mathematical proof of this fact to this date, it 
is expected, on the basis of experiments and numerical simulations that, already at 
moderate Reynolds numbers, the viscous boundary layer detaches from the bound- 
ary and that vortices form in the wake past the immersed body (a phenomenon 
known as a "von Karman street"). Yet there exist perfectly smooth solutions of 
the Euler equations corresponding to the potential fiow of an incompressible fluid 
past a sphere in space dimension 3 — see for instance §10, Problem 2 in [20]. 

Not much is known on Prandtl's boundary layer analysis for this problem from 
the mathematical viewpoint, apart from some positive results based on the Cauchy- 
Kovalevski theorem [10, 11], as well as negative results concerning the Prandtl 
boundary layer equations [14, 17]. There is however the following very interesting 
criterion due to Kato [19]. 

Theorem 1.1 (Kato). Assume that 

[ |u^"(a;)|2rfa; < oo 
Jn 

and denote 

dn^ :={x€n \ dist(a;, dft) < e} . 

Then 




as e — >■ 0. 

In other words, the convergence of the Navier-Stokes solutions to u as e — )■ 
in quadratic mean everywhere in Cl is equivalent to the vanishing of the viscous 
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energy dissipation in a thin layer near the boundary. The convergence of the Navier- 
Stokes solutions to a solution of the Euler equations is therefore a strongly nonlocal 
phenomenon. Notice that, while the Prandtl viscous boundary layer has thickness 
y/e, Kato's criterion involves the vanishing of viscous dissipation in a much thinner 
sublayer, of thickness e. On the other hand, while Prandtl's theory is based on the 
construction of a multiscale asymptotic expansion for Ue, Kato's result is based on 
an energy estimate - so that, in theory, it might happen that Prandtl's construction 
breaks down while converges to u in quadratic mean. (We are however unaware 
of examples of such flows.) 



2. The Inviscid Limit of the Navier-Stokes Equations 
WITH Slip Boundary Condition 

To confirm the role of the Dirichlct boundary condition as a source of difficulties 
in the inviscid limit, we supplement the Navier-Stokes equations with a more general 
class of boundary conditions, known as "slip boundary conditions" , which take the 
form 

(3) Ue-n|gf^=0, e(S(ue) •n)x + Auel^Q = 0, 

where A is a scalar (the slip coefficient), and where E(ue) denotes the deformation 
tensor. In other words, for each vector field v defined on a neighborhood of J7, 

(4) S(t;) := V^v + {V ^vf , 

while the subscript r denotes the tangential component on the boundary: 

(5) Vr ■= V — {v ■ n)n . 
Since • n|g^ = 0, one has 



(i;(-tte) • n)i 



' dUe 



(Vn) • Ut 



so that (S(«,) • n)r\g^ and (^)^ 



differ by a lower order term — of order 0, in 

an 

the sense of differential operators — involving the Weingarten endomorphism Vn 
acting on the tangent space of dO,. 

Henceforth we assume that the slip coefficient A depends on the Reynolds number 
e~^, and denote it by A = Ae. Thus 

a) if Ae > Aq > for all e > 0, one expects that the slip boundary condition (3) 
should be asymptotically equivalent to the Dirichlet condition, thereby leading to 
the same difficulties as regards the inviscid limit; 

b) if Ae = 0, the slip boundary condition (3) reduces to the Navier full slip 
condition 

(6) Me -nj^j^ = 0, (S(?ie) • n)^|gj^ = 0. 

As is well known, the Navier full slip condition prevents the detachment of viscous 
boundary layers, so that the inviscid limit of the Navier-Stokes equations supple- 
mented with that condition is indeed governed by the Euler equations, whenever 
there exists a smooth solution of these equations — see for instance [3] or step 2 in 
the proof of Theorem 4.1 in [21]. 
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In fact, a more general result is true. Consider the Navicr-Stokes equations 
supplemented with the slip boundary condition 

diVx Ue = , 

dtUe + diVx{Ue (X) Me) + VxPe = ^^xUe , O; S O , t > , 

"^■"^100 = 0' e(S(uJ-n)^ + AeUe|aQ = 0, 



(7) 



lt=0 

The existence of global weak solutions "a la Leray" of (7) is classical: see for instance 
Theorem 2 in [18]. 

Henceforth wc use the classical notation T-L{Vl) to designate the set of solcnoidal 
vector fields on O that arc tangent to d^l and have finite square mean in i.e. 

niVl) ■= {v e L^{9.) I div-y = and i; • nj^^^ = 0} . 

Theorem 2.1 (Bardos-Golsc-Paillard [7]). Let S Cl{^) and for each e > 
let Ue be a Leray solution of (7). Assume that the Euler equations (2) have a local 
classical solution u G C^([0, T] x Cl) defined for some T > and satisfying 

(8) / / \u{t,x)\^dS{x)dt+ f f \T.{u){t,x)\^dxdt< 00, 



Jo JdQ Ja Jn 

where dS{x) is the surface element on dfl if N = 3 (resp. the length element if 
N = 2). Then 



— =^ sup / — w| {t,x)dx ^ . 
o<t<T Jn 

The proof will serve as a model in the case of the Boltzmann equation; it is given 

in detail below. 

Proof. Leray solutions of (7) satisfy the energy inequality 
(9) 

i / \u^{t,x)\'^dx + e I I ^\Ti{U)^){s,x)\^dxds + \^ I I \u^{s,x)\^dS{x)ds 
Jn Jo Jn Jo Jdn 

< \ [ |u™(x)|2da;. 
Jn 

On the other hand, for each test vector field w e C^(R+ x O) n C(R+; ^(O)) 

Ue-w{t,x)dx— / u"^{x) ■ w{0,x)dx 
I Jn 

+ e / / ^'S{ue) : H{'w){s,x)dxds + Xe / / ■ w{s,x)dS{x)ds 
Jo Jn Jo Jdn 

= / ■ E{w){s,x)dxds + / / VxW '■ — w)'^^{s,x)dxds 
Jo Jn Jo Jn 

where 

(10) E{w) ■.= dtW + w-VxW- 

Therefore, with the identity 
d 



dt 



^{w^t, x)f dx = / w ■ E{w){t,x)dx , 
Jn 
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we conclude that 

' \u'''{x)-w{0,x)\'^dx 



+e 

(11) 



5 / \Ue - W\'^{t,x)dx - ^ I 

Jq Jn 

n^\j:iue)is,x)\^dxds + Xe I I \u^{s,x)\'^dS{x)ds 
A Jq Jon 

<e / / ^'S{ue) : Ti{w){s,x)dxds + Xe / / ■ w{s,x)dS{x)ds 
Jo Jq Jo Jan 



■ E{w){s,x)dxds — / / ■ {ue — w)^'^{s,x)dxds . 

Jn Jo Jn 

By a straightforward density argument, we may replace w with the solution u of 
(2): thus 

■ E{u){s, x)dxds = — / We • ^xP{s, x)dxds 
_! Jo Jn 

= — / pue{s, x) ■ nxdS{x)ds = , 
Jo Jdn 

so that 
(12) 



; / / \\L{u^){s,x)\'^dxds + K f 
Jo Jn Jo 



\u^ — u\ {t,x)dx+e / / i|I](ue)(s, x)| dxds + Xe / / |Me(s,a;)| dS{x)ds 



an 



< e 



n^S(ue) : S(u)(s, x)dxds + X^ / u^- u{s, x)dS{x)ds 
I Jo Jdn 



+a{u) / / — u\^{s,x)dxds , 
Jo Jn 



where 

a{u) := sup \Vxu{t,x)\. 



0<t<T 



Define 



Qe{t):=e / ^\T.{u,){t,x)\\^{u){t,x)\dx + Xe / \ue{t, x)\\u{t, x)\dS{x) . 
Jn Jdn 

By the Leray energy inequality 

e[ [ ^\i:{ue){s,x)\^dxds + Xe [ [ \u^{s,x)\'^dS{x)ds = 0{1) 
Jo Jn Jo Jdn 

so that, by (8) and the Cauchy-Schwarz inequality, 

r Q,{t)dt = 0{V~e) + 0{^/K)■ 
Jo 

Applying Gronwall's lemma shows that 

/ ^\ue-u\^{t,x)dx Ke'^"^''^ [ Qe{s)ds^O 
Jn Jo 

as e 0. □ 
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3. The Incompressible Euler Limit of the Boltzmann Equation 

3.1. Background on the Boltzmann Equation. The Bohzmann equation gov- 
erns the distribution function of molecules in a monatomic gas, denoted here by 
F = F{t,x,v). We recall that F{t,x,v) is the density with respect to the phase 
space measure dxdv of gas molecules located at the position x € with velocity 
u e at time t. It takes the form 

{dt+v-Va,)F = B{F,F) 

where B{F, F) denotes the Boltzmann collision integral. In the case where gas 
molecules behave like hard spheres subject to elastic, binary collisions, the collision 
integral takes the (dimensionless) form 

B{F,F):= ff {F'Fi-FF,)\{v-v,)-Lo\dv,dio. 

Here we have used the notation 

r F = Fit,x,v), F, = Fit,x,iu), 
^F' = F{t,x,v') , Fl = F{t,x,vl) , 

where v' and v'^ are defined in terms of w, G R'^ and w € by the formulas 

{v' = v'{v,Vt,Uj) := V — {v — V^,) ■ LOLO , 
v'^ = v'^{v,v^,uj) := Vt,+ {v - u*) • WW . 

The Boltzmann collision integral enjoys the following well known properties. 
First it satisfies the local conservation laws of mass, momentum and energy, in 
the following form: if F is rapidly decaying as \v\ oo, then 

(13) / B{F,F)\ Vk\dv = Q, A; = l,2,3. 

Another fundamental property of the collision integral is Boltzmann's H The- 
orem, which can be stated as follows: if F > is rapidly decaying in v while 
InF = 0(|w|") for some n > 0, then 



(14) / B{F,F)\a.Fdv<Q, 

while 

/ B{F,F)lnFdv = Q ^ B{F,F) = Q ^ F = M(p,u,0) , 

JR3 

where the notation M(^p^u,e) designates the Maxwellian distribution with density 
p > 0, bulk velocity u e R^ and temperature ^ > 0, i.e. 

We shall formulate below the corresponding differential relations for solutions of 
the Boltzmann equation. 
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3.2. Boundary Conditions for the Boltzmann Equation. A general class of 
boundary conditions for the Boltzmann equation is of the form 

(15) F(t,x,v)\v-nx\ = / F{t,x,v')v' ■nxlC{x,TZxV,dv') , {x,v)gT_, 

J v'-n^ >0 

where we recall that the notation Ux designates the unit outward normal at the 
point X e dO. while TZx designates the specular reflection: 

TlxV :=v-2v- nxUx , 

and where we have used the notation 

r+ := {(x, v) ednx'R.^\vnx>0}, 

r_ := {{x, v) e dnxll^\v Tlx <o} . 

The measure- valued, scattering kernel IC{x,v,dv') > satisfies the assumptions 
(see §8.2 in [12], especially on pp. 230-231): 
(i) for each x G dfl, one has 



/ 



/C(x, V, dv')dv = dv' ; 
(ii) for each x ^ dO. and each $ € Cc(R^ x R^) one has 



^{v,v'){v ■ nx)+{v' ■ nx)+M(ififl){v)M(ififi){v')dvK{x,'RxV,dv') 

$(-«', -w)(w • nx)+{v' ■ nx)+M(ifl,g){v)M(ifl,g){v')dvJC{x,TlxV,dv') 

R3xR3 

{liK{x,v,dv') is of the form lC{x,v,v') = K{x,v,v')dv' , property (ii) is equivalent 
to the identity 

K{x, V, v') = K{x, v' , v) for a.e. v, v' s.t. v ■ Ux > and v' ■ nx > 

for each x G dfl.) 

(iii) for each x G dfl, one has 

M(^ififi){v)v-nx = / M(ififi){v')v'-nx]C{x,v,dv') , v-nx>0, 

Jv'-n^>a 

if and only if 6 = 9^ {x) , where (x) is the temperature of the boundary dfl at the 
point X. 

With these notations, the case of diffuse reflection corresponds to 

M(ifl^0^){v)\v ■ UxW 



}C{x,v,dv') := 



/ M(^ifi^0^){u)\u- nx\du 

Ju-nx<0 



The case of an accommodation boundary condition (i.e. the so-called Maxwell-type 
condition considered in [7]) corresponds to 

■A^(i,o,e„)(w)|w • nx\dv' 



K{x, V, dv') := (1 — a{x))6{v' — v) + a{x) 



/ M(^ifi^e^){u)\u-nx\du 

Ju-nx<0 



with < a(.T) < 1. The case a = 1 corresponds to diffuse reflection, while the case 
a = corresponds to specular reflection — notice that specular reflection does not 
satisfy assumption (ii). 
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The class of reflection kernels described here also includes the Cercignani-Lampis 
model, as well as all Nocilla models. We refer the interested reader to section 8.4 
in [12] (on pp. 235-239) for more information on this issue. 

However, this is by no means the most general class of admissible boundary 
conditions for the Boltzmann equation. Indeed, the boundary condition (15) is 
local in the position variable x. There are also models of gas-surface interaction 
involving boundary conditions satisfied by the distribution function of the gas that 
are nonlocal is both the position and velocity variables: see for instance [23, 24]. 
This new class of boundary condition is not covered in the present study. 

3.3. The Incompressible Euler Scaling for the Boltzmann Equation. The 

incompressible Euler limit of the Boltzmann equation is based on three different 
scaling prescriptions (see for instance [4, 5]). A careful description of the dimen- 
sionlcss form of the Boltzmann equation can be found in section 1 of [6], as well as 
in sections 1.9 and 1.10 of [30]. 

First, as in all fluid dynamic limits of kinetic models, one assumes a strongly 
coUisional regime. In other words, the collision integral is scaled as 

B{F, F) = -^B{F, F) , with e < 1 and g > . 

(See [6] or [30] for the physical meaning of e and q.) 

Next, the incompressible Euler limit holds on a long time scale, leading to the 
introduction of a new time variable t — t/e: in other words, we set 

F{t, X, v) = F^{i, .X, v) = F^{et, x, v) . 

With these assumptions, the scaled Boltzmann equation becomes 

{edi + vS/^)F, = -^^B{F,,F,). 

A third and last scaling assumption used in the incompressible Euler limit of 
the Boltzmann equation is that the corresponding gas flow is kept in a low Mach 
number regime. Specifically, we assume that the Mach number is of order IVb, = e, 
and this is done by seeking the distribution function F^ in the form 

Fe = A^(i,o,e„)Ge with = 1 + eg^ 

(where it is implicitly assumed that ge — 0(1))- 

We draw the reader's attention to the fact that the parameter e considered here 
has a different meaning than in the inviscid limit of the Navier- Stokes equations. 
Indeed, in the situation considered here, we shall see that the reciprocal Reynolds 
number is Re""'^ = e'^ — instead of e as in the previous section. 

For simplicity, we henceforth drop all hats in the scaled Boltzmann equation and 
consider the initial boundary value problem 
(16) 

' {edt + vV,)F, = -^BiF„F,), {x,v)€nxR\ 
' F^(t,x,v)\v-nx\ = / Fe{t,x,v')v'-nxJCe{x,TZxV,dv') , {x,v)er_, 

Jv'-nx>0 
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assuming that the boundary temperature 6'™ = 1. The reflection kernel /Ce satisfies 
the properties (i)-(ii) above and possibly depends on the scaling parameter e. 
Henceforth we denote 

M:=7W(i,o,i). 

Let us return to the conservation laws of mass, momentum and energy satisfied 
by the Boltzmann collision integral. If is a smooth solution of (16) rapidly 
decaying as \v\ — >■ +cxd, then it satisfies the system of differential identities 



(17) 



edt / F^dv + diva: / 
iR3 Jb 



vF^dv = , 



vF^dv + divx 



R3 



R3 

V ® vF^dv = , 



R3 



edt / ^\vf F^dv + div a: / v^\vf 

^ JR3 JR3 



■F,dv = . 



These identities can be viewed again as the differential form of the local conservation 
laws of mass momentum and energy that are classical in continuum mechanics. 

3.4. Main result. Our main result in this paper is an analogue of Theorem 2.1 

for the incompressible Euler limit of the Boltzmann equation. 

Theorem 3.1. Let m*" S C^(0;R^) satisfy diVxU^"' = and u ■ = 0, and 
assume that the Euler equations (2) have a local classical solution u G ([0, T] x fi) 
defined for some T > and satisfying (8). For each e > 0, let he a solution of 
(16). Assume that the reflection kernel K.^ satisfies, in addition to the properties 
(i)-(iii) listed above, 



(18) 



v'^dv' — I VrlCe{x,v,dv')dv 

Jv-nx>0 



< ae{x)\v'\dv' 



and 

(19) }C^{x,v,dv') > ae{x)P^{x)V2^M{v){v ■ n^)+dv' 

for each x G dfl and a.e. v, v' G such that v ■ Ux > and v' ■ nx> 0, where 



(20) 

and 
(21) 



< oif^{x) < 1 and 



- sup ae(a;) 

e xG3n 
\m\<R 



as € 







< inf inf Pe{x) < P,{x) < 1 

|a!|<R 



for each R> 0. Then, for each T > and R> 

vF^it, X, v)dv — u(t, x) 



lxl<R 



R3 



dxdt — >■ as e — >■ . 



This result was stated and proved in [7] in the particular case where the reflection 
kernel K^^ corresponds with Maxwell's accomodation condition at the boundary with 
accomodation parameter a^. 

In the case where the Boltzmann equation is set on a domain without boundary 
— i.e. if a; runs through the Euclidian space R"^ or the periodic box T'^ — the 
first rigorous derivation of the incompressible Euler equations from the Boltzmann 
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equation was obtained by L. Saint-Raymond [28] following the relative entropy 
method sketched in [9, 22]. (For the simpler case of the BGK model, sec also [27].) 

The proof of Theorem 3.1 is a slight generalization of the one in [7] and is based 
on the same relative entropy method as in the work of L. Saint-Raymond [28]. It 
can also be viewed as the analogue for the Boltzmann equation of the proof of the 
inviscid limit Theorem 2.1. 

A final word of caution is in order. The statement of the theorem is left delib- 
erately vague about the notion of solution of the Boltzmann equation (16) to be 
considered. This is a rather technical matter, to be discussed later. 

3.5. The relative entropy inequality. The first important property of the Boltz- 
mann equation used in the proof is the following variant of Boltzmann's H theorem 
that can be viewed as an analogue of the Leray energy inequality (9): 

\H{F,\M){t) + ^ f I r,{s,x)dxds 
(22) ' ' 

+ 4 / / 'DGe{s,x)dxds < \HiF,\M){0) 

where the notation H{Fi\F2) designates the relative entropy defined as follows. 
Let Fi > and F2 > a.e. on il x designate two measurable functions, then 

H{Fi\F2) ■■ = JJ^ (^1 In (|r) - ^1 + ^2) {x,v)dxdv 

[Fi \ 
h\ — — 1 1 ^2(0;, v)dxdv , 

OxR3 \F2 J 

where h is the function defined as follows: 

h: [-1, 00) 9 z 1-^ h{z) := (1 + z) ln(l + z) - z €'R+ . 

(Notice that, since the integrand h — 1^ F2 is a nonnegative measurable func- 
tion, the relative entropy H[Fx\F2) is always a well-defined element of [0, -|-cxd].) 

The two other quantities in (22) are the entropy production rate per unit volume 
Vc and the Darrozes-Guiraud information at the boundary DQ^, whose definition 
is recalled below. 

The entropy production rate per unit volume for the Boltzmann equation is 
Ve-- [ B{F„F,)\nF,dv 
and can be put in the form 

Ve= III r(^^-l\F,F,,\{v-v,)-uj\dvdv,du:, 

JJjR3xR3xS2 \i'ei'e* J 

where r is the following function: 

r : (—1, 00) 9 z H- > r{z) :— z ln(l + z) G R+ . 
As for the Darrozes-Guiraud information, it is defined as 

Vg^{t,x):= j h(^-l^{t,x,v)v-n^Mdv, x€dn,t>0. 
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Observe that, for each x E dft and each v G R'* such that v ■ Ux > 0, 

e . , /^ jv' ■nx)+Miv')IC,ix,v,dv') 
{v.n,)+M{v) 

is a probability measure on (by property (iii) of /C^), and that 

^{t,x,nxv) = /nl^y, ^/(^'^''") ' for each {x,v) e r+ and t > . 



Since h is convex, it fohows from Jensen's inequahty that 
On the other hand, by property (i) of the reflection kernel /C^, one has 



Hence the Darrozes-Guiraud information satisfies 

VgS,x) = J^^ l^h (^^-?j-h (^^//,,,, § - l))) M{v ■ n.)+dv > 0. 

(See for instance Theorem 8.5.1 on p. 240 in [12].) 

As a consciqncncc of the local conservation of mass, i.e the first equality in (17) 
one finds that, for each scalar test function (p e C^(R+ x Q) and each t > 0, 

ell Fc{t,x,v)4>{t,x)dxdv — e // Fe{0,x,v)(j){0,x)dxdv 



(23) + / // F^(s, x,v)(j){s, x)v ■ nxdS{x)dvds 

Jo JJdQxn^ 



= // F^{s,x,v){edt + V ■ 'Vx)(t>is,x)dsdxdv . 
Jo JJnxR3 



an • 



Let now w G Cl{[0,T] x fl;R3) satisfy 

(24) diVa; w = , and w ■ n\ 

As a consequence of the local conservation of momentum, i.e. the second identity 
in (17), one finds that, for each t > 0, 

Fe{t,x,v)v ■ w{t,x)dxdv — e // F^{0,x,v)v ■ w{0,x)dxdv 

(25) + / // F^{s,x,v){vr ■ Wr)v ■ nxdS{x)dvds 

Jo JJdnxR^ 

= // F^{s,x,v){edt +v ■'Vx)iv ■ w){s,x)dsdxdv . 

Jo iioxR!* 

Finally, observe that 



(26) 



^^(i"e|7W(i,.».i)) = H{F,\M) + [[ FJn ( ) dxdv 

= H{F,\M) + [[ F, (ie2|u;|2 -ewv) dxdv. 
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Putting together (26), (22), (23) and (25), we find that, for each test velocity 
field w e Cc([0,T] x OiR^) satisfying (24) and each t>0 

^H{F,\M^i.,^.X)m < ^//(Ff |M(i,,^(o,),i)) 

1 /■ ^ f'^ f 

P^{s,x)dxds ^ / / VQ^{s,x)dS{x)ds 

e Jo Jas 



Jn t Jo Jdn 



1 



C27") 2 / // (i) — ew(s, x))*' : Vxw(s, a-')^£(*"j ^'j ^O'^sjdwds 

— - / // {v — ew{s,x)) ■ E{w){s,x)Ff{s,x,v)dvdxds 
e Jo J JnxRW 

~2 f ff F^{s,x,v){vr ■ Wr){s,x){v ■ nx)dS{x)dvds , 
e Jo JJanxRW 

where i?(w) has been defined in (10). 

This inequality is the analogue for the scaled Boltzmann equation of the in- 
equality (11) used in the proof of Theorem 2.1. Indeed, for each pair ui,U2 of 
(measurable) vector fields in ft such that 

/ \ui{x)\'^dx+ / \u2{x)f dx < +00 , 
Jn Jn 

one has 

^-ff(A^(l,e«2,l)l-^(l,e«i,l)) 5 / \UI-U2\'^{x)dx. 

e Jn 
Therefore, the scaled relative entropy 

-^H{F^\A4(^i^ew,i)) is the analogue of ^ / \ue — u\'^{t,x)dx , 
e ' ' Jq 

the entropy production 

1 r* r 

Veis, x)dxds 



_L/* 

is the analogue of the viscous dissipation 



( j \\L{u^){s,x)\'^dxds, 
Jn 



— with the notation and scaling assumption used in sections 1-2 — while the 
Darrozes-Guiraud information 

4 / / Vg^{s,x)dS{x)ds 
e Jo Jan 

is the analogue of the boundary friction 



'0 Jan 

appearing in section 2. Likewise the term 



Ae / / \ue{s,x)\'^dS{x)ds 
Jo Jon 

ewise the term 

-2 /" 11 {v — e'w{s,x))^'^ -.V xw{s,x)F^{s,x,v)dxdvds 
e Jo JJoxR'v 

is the analogue of 

/ / y xW : {u^ — w)®'^{s,x)dxds , 
Jo Jo, 
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and we seek to control it in terms of the scaled relative entropy 

^iJ(i^,|7W(i,,^,i)) 

in order to conclude with Gronwall's inequality. This is precisely what is done in 
[28], and we shall not repeat this (difficult) analysis here. 
What remains to be done is to control the boundary term 

1 ''' 

c2 



F^{s, X, v){vr ■ Wr ){s, x){v ■ nx)dS{x)dvds 

JJanxTL^ 



of indefinite sign by the Darrozcs-Guiraud inequality, possibly up to some asymp- 
totically negligible quantity. This step can be viewed as the analogue for the Boltz- 
mann equation of the Cauchy-Schwarz inequality used to control the boundary 
term 

Ac / / ■ u{s,x)dS{x)ds 
Jo JdQ 

by the boundary friction 

Xe I I \u,\^{s,x)dS{x)ds. 
Jo Jon 

3.6. Proof of Theorem 3.1. The key argument in the proof of Theorem 3.1 is 
given by the following proposition. 

Proposition 3.2. Under the same assumptions as in Theorem 3.1, for each vector 
field w G C^(R+ X f2;R^) such that dwxW = and w ■ n\Q^ = 0, satisfying 
supp(w;) C R+ X K, where K is a compact subset ofTL^, one has 

/ // F^{s,x,v){vr ■ Wr){s,x){v ■ nx)dS{x)dvds 
Jo JJonxn^ 



< h^ + ^':?:;^^ ^ / / ^^GM^)dxds 



Ne h{l/2) J Jo JoanK Pe{x) 

+ 2C{w,Ny [ [ ^^iM /" F^{s,x,v){vnxfdvdxds 
Jo JdnnK £ Jr^ 

for each e S (0, 1), each N > 1 and each t > 0, where 

(28) C{w,N):=^ [ (e^ll'"ll^~l''l-7V||w||Loo|t;|-l)(t;-n^)+Mrft;. 
Set 

(29) F, = MG, , G, = l + eg,. 

Using the boundary condition in (16) and the substitution v TZxV fov v ■ < 
shows that 

/ // F^{s,x,v){vr ■ Wr){s,x){v ■ nx)dS{x)dvds 
Jo JJanxRs 

= / // {Fe{s,x,v) — Fc{s,x,TZxv){vt ■ Wt){s,x){v ■ nx)+dS{x)dvds 
Jo JJanxii^ 



14 



F. GOLSE 



notice that TZxVt = Vt- Thus 

ft 



III 

Jo J Jo 



anxR3 



F^{s,x,v){vj- ■ Wt){s,x){v ■ nx)dS{x)dvds 



{ge{s,x,v) — ge{s,x,TZx'i>){vr ■ Wr){s,x){v ■ Tlx) +M {v)dS {x)dvds 

JJcKIxh^ 

{9e - {lJ'l,v,9e)){s,x,v){vr ■ Wr){s,x){v ■ nx)+M{v)dS{x)dvds . 

lO JJdQxR^ 

Since ^ is a probability measure acting on the velocity variable only, 

9e - {lJ'l,v,9e) =9e- A^ge - (mI,^, 5e - A^^e) , 

where 

(30) A,<A : 

Therefore 



= e 



:= V2Tr / (}>{v){v ■ nx)+M{v)dv . 

J J Fe{s,x,v){vT ■ Wt{s,x,v)){v ■ nx)dS{x)dvds 

{9<L - Kge - {lA,v^9e - Axge))ivT ' Wr)iv ■ nx)+M{v)dv 

I // ae{x)\ge- Axge\{s,x,v)\vr\\wr\is,x){v ■nx)+M{v)dv 
J Jr+ 



< e 



by assumption (18). 

At this point, the proof of Proposition 3.2 is done in two steps. 
The first step is summarized in the following lemma. 

Lemma 3.3. Under the same assumptions as in Proposition 3.2 and with the 
notation (29), for each {t,x) € (R+ x 9i7)nsupp(w), each N > 1 and cached (0,1), 

^^r||wr|(^^ • nx)+M{v)dv 

Jr3 

< [ {h{eg,)-h{eAx{gemv-nx)+M{v)dv + NC{w,N)e^Ax{F,), 

where A^ has been defined in (30) and C{w,N) in (28). 

Proof. The proof follows that of Lemma 4.6 in [7]. Pick zq > —1, and define, for 
each z > —1 

l{z ~ zq) := h{z) - h{zo) - h'{zQ){z - zq) . 
We recall that the Legendre dual of h and I are 

h*{p) = eP-p-l and I* (p) = (1 + za)h* (p) = (1 + zo)(eP - p - 1) . 

Applying Young's inequality (see [6] , especially section 3 there) to I with z = egc 
and zo = eA^ige) 

Ne'^lge - Ax{ge)\\vr\\wr\ < l{e{ge - Ax{ge))) + l*{Nesign{g^ - Ax{ge))\vr\\wr\) 
= l{e{g, - AM)) + (1 + eA,(ff,))(e^^l'"ll''l - Ne\w\\v\ - 1) 
< l{e{g, - AM)) + e2(l + eA,(ff,))(e^l-ll''l - N\w\\v\ - 1) , 
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whenever < e < 1. (Indeed, for each z > and e G (0, 1), one has 



n! 



n>2 

Z 

n\ 



n>1 



see also (3.14) in [6].) 
Therefore 

\ \9(.- l^x{g€)\{s,x,v)\vr\\wr\{s,x){v ■nx)+M{v)dv 

(31) 

< \ l{e{g, - AM)){v ■ n,)+M{v)dv + NC{w, N)e''K{F,) . 
Moreover, since is constant in v and the average under a probabihty measure, 
/ h'{€AM){9e - Ax(5e))(v • n^)+Mdv 

>/r3 

= h'{eA:,{ge)) / {ge - Ax{ge)){v ■ n:,)+Mdv 

JR3 

= h\eAM)^^x{ge - AM) = , 

so that 

l{e{ge - A^{ge))){v ■ n^)+M{v)dv 



I' 

JR3 



= / {h{ege)-h{eAM)){vn^)+M{v)dv. 

iR3 

Substituting this last integral in the right hand side of (31) leads to the inequality 
in the statement of Lemma 3.3. □ 

It remains to relate the integral on the right hand side of the inequality in Lemma 
3.3 to the Darrozes-Guiraud information. 

Lemma 3.4. Under the same assumptions as in Proposition 3.2 and with the 
notation (29), for each {t,x) G R+ x dQ and each e G (0, 1), one has 

a,{x)l3,{x) [ {h{eg,) - h{eAM)){t, x, v){v ■ n^)+M(y)dv < Vg,{t, x) . 

JR3 

Proof. By (19), the reflection kernel K, can be put in the form 

JC{x, V, dv') — ae{x)l3e{x)V2Tr{v ■ nx)-f-M{v)dv + (1 — ae{x) l3e{x)) Ce{x , v, dv') . 

Since < ae(a;)/3e(x) < 1, the kernel £e{x,v,dv') so defined satisfies the properties 

(i), (ii) and 

(iii') for each x G dSl, 

M{v)v-nx= / M{v')v' ■nxCe{x,v,dv') , v-nx>Q. 

Jv'-n^>0 

Define i'x,v by analogy with iJx,v' 

{v' ■ nx)+M{v')Ce{x, V, dv') 



{v ■ nx)+M{v) 
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SO that the boundary condition satisfied by the solution Fe of (16) takes the form 

^{t,X,n^v) = aeix)l3e{x)A^ ^ + (1 - ae{x)Pe{x))(^iyl^, ^^i*^^^^) ' 

for each (a;, v) G r+ and t>0, 

or equivalently 

ge{t,x,nrov) = ae{x)Pe{x)A^{g^){t,x,v) + (1 - ae{x)l3^{x)){vl .^,g^{t,x, •)) , 

for each {x, v) € r_|_ and t> . 

Thus 

'DQ^{t,x) = J h{eg^){t,x,v)v ■ nxM{v)dv 

= J {H<^9€){t,x,v) - h{eg^){t,x,nxv)){v ■ nx)+M{v)dv 
which can be put in the form 

= -^^-i^i.^9e) - h{a,/3,{x)AM + (1 - (^ePe{x)){l^l^,,ge))) 

By convexity of h 

h{aePe{x)A^{g,) + (1 - a,l3,{x)){iy;^ ,^, g,)) 

< a,p,{x)h{K{g,)) + (1 - a,P,{x))h{{vl^^,g,)) 

so that 

Vg, > -;^a,p,{x)K^{h{eg,) - h{AM)) 

+ vfe(l - a,0,ix))A,ihieg,) - h{{ul^,g,))) 
Since the kernel Ce satisfies property (i), one has 

A,{h{eg,))=A,{{ul^^,h{eg,))) 

so that 

Ax(/i(effe) - h{{ul,,g,))) = A^{{yl^,h{eg,) - h{{ul„g,)))) > 0, 

where the last inequality follows from Jensen's inequality applied to the probability 
measure ^ and the convex function h. 

Since < ae(a;)/3g(a;) < 1 for all x G dSl, we conclude that 

m > ^a,^,{x)A^{h{eg,) - h{A^{g,))) 

which is precisely the sought inequality. □ 

Finally, we control the outgoing mass flux at the boundary exactly as explained 
in Lemma 4.7 of [7]. 

Lemma 3.5. Under the same assumptions and with the same notation as in Propo- 
sition 3.2, for each {t,x) G x dfl and each e,ri G (0, 1), one has 

a^ix) [ F,{t,x,v){v ■ n^)+M{v)dv < \ , Vg,{t,x) 

JR3 



27r(l - ri) Jr3 
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Proof. Following the argument in Lemma 4.7 of [7] shows that 

Fe{t,x,v){v ■ nx)+M{v)dv 



R3 



^TT-T / (GJn(^-%—]-G, + A,{G,)]{t,x,v){vn,)+M{v)dv 



+^^=7! Fe{t,x,v){v-na;fdv 

\/27r(l - r?) Jr3 

(/i(ege) - h{eAx{g^))){t,x,v){v ■ nx)+M{y)dv 
V27r(l - r/) 



We conclude by applying Lemma 3.4 to the first integral on the right hand side of 
the last equality above. □ 

Putting together the inequalities in Lemmas 3.3, 3.4 and 3.5 with V = ^ leads 
to the estimate in Proposition 3.2. 

Next we return to (27), and observe that the boundary term can be absorbed 
by the Darrozes-Guiraud information, as follows. 

Pick a compactly supported vector field w of class such that div^; w = 0, 
w ■ n\g^ =0 and let be a compact subset of R'^ such that supp(t(;) c R+ x K. 
Set 

CM. { X) 

:= sup — and /?* := inf I3e{x) > 0; 

— see assumptions (20)-(21). Pick A'' > 1 such that iV/3, > 1. It follows from 
Proposition 3.2 that 

\H{F,\M^,,^,l){t) < 4i?(^^ri^^l,e^(0,.),l) - ^ / / 'Pe{s,x)dxds 

£ £ ^ ^ Jo Jn 

H1/2)P. ' ) I l^^0.is,x)dSix)ds 

(32) ~~2 [ II {v - ew{s,x))®'^ ■.Wxw{s,x)F^{s,x,v)dxdvds 

— - / // {v — ew{s,x)) ■ E{w){s,x)F^{s,x,v)dvdxds 
e Jo J JQxn^ 

ft r r 

^2, 



+2C(w,N)re // F^{s,x,v){v ■nxydS{x)dvds. 

Jo JJ(aanA:)xR3 

At this point, we recall that solutions of the initial boundary value problem (16) 
Boltzmann equation satisfy the bound 

sup / / / F^{s,x,v){v ■ nx)'^dS{x)dvds < 00 

00 Jo J J{dnnK)x'R.^ 

for each t > and each compact K (ZR^ — see [25] and the proof of Theorem 4.1 
d) in [7]. 
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Since — J- as e — )■ by our assumption (18)-(20) on the reflection kernel /Ce, 
and 

\ /i(l/2)/3, J ^ 

for all small enough e > 0, the inequality (32) is of the same form as the inequality 
stated as Theorem 5 in [28]. One then concludes by the same argument as in [28]. 

4. Conclusion and Final Remarks 

In this work, we have derived the incompressible Euler equations from the Boltz- 

mann equation in the case of the initial boundary value problem, for a class of 
boundary conditions at the kinetic level of description that is more general than the 
case of Maxwell accomodation considered in [7]. However, the asymptotic regime 
for these boundary conditions considered here is the same as the one considered in 
[7]: condition (18)-(20) means that, in the vanishing e limit, the tangential momen- 
tum of gas molecules tends to be conserved. In the case of Maxwell's accomodation 
condition at the boundary, this is the case precisely when the accomodation pa- 
rameter vanishes, so that Maxwell's condition approaches specular reflection. The 
(unrealistic) specular reflection condition is the analogue for the Boltzmann equa- 
tion of the (equally unrealistic) Navier full slip boundary condition, for which the 
inviscid limit of the Navier-Stokes equation is known to be described by the Euler 
equation. In that sense, Theorem 3.1 is the analogue for the Boltzmann equation 
of the (much simpler) Theorem 2.1 already established in [7]. 

There is no single notion of accomodation parameter for boundary conditions 
defined by a reflection kernel as in (15), as in the case of Maxwell's condition. 
There is however a notion of accomodation parameter that can be defined for each 
mechanical quantity attached to a gas molecule impinging on the material surface; 
see the discussion in section 8.3 in [12], and especially formulas (3.4) and (3.11) 
there. The interested reader is invited to compare that definition of accomodation 
parameter with the condition (18)-(20) used in Theorem 3.1. 

We have already stressed the striking analogy between our proofs of Theorems 
2.1 and 3.1. In order to go further in the analysis of the incompressible Euler limit 
of the Boltzmann equation in the presence of material boundaries, it seems natural 
to investigate the similarities between the various terms that appear in the energy 
inequality (11) in the Navier-Stokes case, and the relative entropy inequality (27) 
in the Boltzmann case. The entropy production rate e~^~'^Ve in (27) is clearly the 
analogue of the viscous dissipation rate e|S(U(;)p (with the notation and scaling 
assumptions used in sections 1 and 2: sec Proposition 4.6 in [6]. and especially 
formula (4.18) thcrc^. The analogy with the Darrozos-Guiraud information e^^DQe 
is already somewhat less clear. The boundary friction term Ag|Mgp appearing in 
(11) is partly analogous to the Darrozes-Guiraud but only partly so, since the 
Darrozes-Guiraud information also controls the departure of the distribution func- 
tion of gas molecules impinging on the material boundary from thermodynamic 
equilibrium. For instance, the Darrozes-Guiraud information term is present ex- 
cept in the case of specular reflection of gas molecules at the boundary, whereas 

^This formula is established under the assumption of the Navier-Sokes scaling corresponding 
with g = in the notation of the present paper; extending it to the Euler scaling is straightforward. 
The reader is also invited to pay attention to the definition of q in that reference, whose meaning 
is different than in the present paper. 
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the boundary friction term vanishes identicahy if the solutions of the Navier-Stokes 
equations are assumed to satisfy the Dirichlet boundary condition. 

If the family of solutions of the Navier-Stokes equations (1) or (7) converges 
to the solution of the Euler equations (2) as e — >■ 0, then 

e|S(we)r ^ in L'^iiO,T) x O) 

where T is the life-time of the Euler solution, and, in the case where the Navier- 
Stokes solutions satisfy a slip boundary condition with slip coefficient Ae, 

Ae|ue|af^|^ in L^((0,T) x dQ) 

as e — > 0. (Indeed, the kinetic energy of the fluid is an invariant of the motion 
under the Euler dynamics.) Likewise, if the family of solutions of (16) satisfies 

If 1 

- / vFedv u and -^H{F^\M(^i^eu,i)) 

as e 0, where u is the solution of (2), then 

-J+^Ve -S^ in L^((0,T) x ft) and ^VQ^ in L^{{0,T) x dSl) . 

In the case where the Navier-Stokes solutions satisfy the Dirichlet condition, it is 
easily seen that the Navier-Stokes solutions converge to the Euler solution if and 
only if e(E(u,) ■ n)r ^ in V'{{0,T) x d^l)). In Kato's paper [19], this term is 
controlled by only the viscous dissipation term, by a very clever argument involv- 
ing the construction of a boundary layer different from the Prandtl construction 
together with Hardy's inequality. Whether there is an analogue of this argument 
for the Boltzmann equation remains an open question at the time of this writing. 

Another issue related to the ones considered here is the incompressible Navier- 
Stokes limit of the Boltzmann equation with boundary condition of the type (15). 
At the formal level, this question is investigated in detail in [30]. However, the dis- 
cussion in [30] leaves aside the special case where the Navier-Stokes limit leads to a 
slip boundary condition. In the case of the BGK model with Maxwell's accomoda- 
tion condition at the boundary, this problem has been treated (at the formal level) 
by Aoki-Inamuro-Onishi in [1], by the Hilbcrt expansion method, completed with 
appropriate Knudsen layer terms. In the case of the Boltzmann equation, under 
a scaling assumption leading to the evolution Stokes equations, the slip boundary 
condition has been obtained by Masmoudi and Saint-Raymond [25] by a rigorous 
argument involving the same moment method as in [5, 16] together with the weak 
(i.e. in the sense of distributions) formulation of the slip boundary condition. The 
Navier-Stokes analogue of this result has been obtained (at the formal level) in 
section 3 of [7]. 

While both approaches lead to the same value of the slip coefficient, a natural 
question is to study the Navier-Stokes limit of the Boltzmann equation with the 
most general boundary condition of the form (15), and to identify under which 
condition(s) on the reflection kernel /C, the limiting velocity field satisfies the slip 
boundary condition. In particular, one should check that the Aoki-Inamuro-Onishi 
theory and the weak formulation of the momentum equations lead to the same 
boundary condition in the vanishing e limit. We hope to return to these questions 
in a future publication [15]. 

Finally, we conclude with remarks of a more technical nature, on the notion 
of solution of the initial boundary value problem for the Boltzmann equation. A 
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careful examination of tfie proof of Theorem 3.1 reveals that, in addition to the 
properties of solutions already used in [28] in the absence of material boundaries, 
all that is needed is the relative variant of Boltzmann's H theorem (22), the weak 
form of the local conservation of mass (continuity equation) (23) and the weak form 
of the local conservation of momentum (25). Recently, Mischler [26] constructed 
global, renormalized (in the sense of DiPerna-Lions [13]) solutions of the Boltzmann 
equation in (16), satisfying (22) and (23). However, the boundary condition in (16) 
should be replaced with the inequality 



Equality in the boundary condition is known at least in the case of Maxwell's ac- 
comodation condition — see Remark 6.4 in [26]. The reason for these differences 
with what would be normally expected from classical solutions of (16) is that Mis- 
chler's solutions are obtained by a very delicate compactness procedure in some 
weak topology. Whether these solutions are uniquely determined by their initial 
data, or even satisfy the weak form of the local momentum conservation (25) re- 
mains yet unknown. On the other hand, obtaining classical solutions of (16) for all 
e > and for all compactly supported and solenoidal initial velocity fields of 
class in the domain fl and tangential on dfl remains an open problem at the 
time of this writing. Since the only properties of solutions used in our proof of the 
Euler limit are (22), (23) and (25), one could hope that some future refinement of 
Mischler's theory could be enough for our purposes. 

On the other hand, the Euler limit stated in Theorem 3.1 can be formulated in 
terms of dissipative solutions of the incompressible Euler equations (see [21] for this 
notion, and [8] for its extension to the initial boundary value problem). The main 
advantage of such a formulation is that a) the definition of dissipative solutions 
corresponds exactly with the limiting form of the relative entropy inequality (27) 
as e — > 0, and b) for any initial velocity field with finite kinetic energy that is 
solenoidal and tangential on dfl, there exists at least one dissipative solution of the 
incompressible Euler equations that is defined for all times — see Proposition 4.2 
on p. 156 in [21]. Besides, should the Euler equations have a classical solution, it is 
known that all dissipative solutions with the same initial data must coincide with 
that classical solution — see Proposition 4.1 on p. 155 in [21]. We have refrained 
from using the notion of dissipative solutions in the present paper for the sake 
of simplicity, and refer the interested reader to [7] where the result analogous to 
Theorem 3.1 in the special case of Maxwell's accomodation condition is formulated 
in terms of dissipative solutions of the incompressible Euler equations. 
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